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Abstract 

We present a continuous time generalization of a random walk with complete memory of its 
history [Phys. Rev. E 70, 045 101 (R) (2004)] and derive exact expressions for the first four moments 
of the distribution of displacement when the number of steps is Poisson distributed. We analyze 
the asymptotic behavior of the normalized third and fourth cumulants and identify new transitions 
in a parameter regime where the random walk exhibits superdiffusion. These transitions, which are 
also present in the discrete time case, arise from the memory of the process and are not reproduced 
by Fokker-Planck approximations to the evolution equation of this random walk. 
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I. INTRODUCTION 



Recently, Schiitz and Trimper obtained the exact moments of the probability distribu- 
tion function (PDF) of total displacement of a discrete time random walk with memory 
of its entire history Q]. They coined the term "elephant random walk" (ERW) for this 
model, referring to the myth that elephants have perfect memories. In the ERW, long-range 
memory effects enter the system dynamics through the microscopic evolution equation as 
a displacement- and time-dependent bias. This microscopic approach to modeling memory 
contrasts sharply with the largely established practices of introducing correlated random 
noise terms in Langevin equations [2 L memory kernels in generalized Fokker-Planck equa- 
tions or combinations of both at a coarse-grained scale. 

In its simplest formulation, the discrete time ERW takes place on a one-dimensional 
lattice with lattice constant A and each step occurs at time intervals 1/A. The memory 
of this process is controlled by the initial bias (3 = 2q — 1 and the correlation parameter 
a — 2p — 1. Here, q is the probability of taking the first step in the positive direction while 
p is the probability of taking succeeding steps in the same direction as a randomly chosen 
displacement in the walker's history. The corresponding microscopic evolution equation for 
the probability Pn+i{u) of the walker having a total displacement y from the origin after 
N + 1 steps is therefore 



Px(y) = ~(1 - mv + A) + ^(1 + 0)6(y - A), 



N+l 



(y) 



+ 



®{y + A) 

NA 

a (y - A) 

NA 



PN(y + A) 

Pn(v - A), N>1, 



(1) 



where 8{y) is the usual Dirac delta distribution. Details of the stochastic rules governing 
the movement of the elephant random walker may be found in 

The sequence of displacements in an ERW can also be pictured as a sequence of binary 
events. Such a perspective would be of great use in the study of two-level systems where 
successive occupation probabilities are linear functions of the fraction of occupants in the 
respective levels. Some examples are two queues with arrivals biased toward the shorter 
one and two competing states or agencies with recruitment favoring the more popular one. 
Such an interpretation has been successful in the study of correlated bit strings observed in 
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reduced sequences of DNA bases in genomes, letters in written text, and losses and gains in 
the stock market j?i 0] • 

An analysis of the scaling properties of the mean and mean squared displacement of the 



PDF of total displacement reveals two distinct transitions in the dynamics of the walker 



1. When the parameter a is negative, the sequence of displacements is anticorrelated so 
that the walker is localized about the origin. Conversely, positive correlations lead to 
a net drift in the direction of the initial bias when a > 0. An uncorrelated random 
walk is recovered when a is exactly zero. 

2. When a is less than 1/2, the mean squared displacement scales linearly with the 
number of steps N as in a normal diffusive process. On the other hand, when a > 1/2, 
the mean squared displacement scales asymptotically as a power law that depends on 
a and results in super diffusion. When a is exactly 1/2, the process is marginally 
superdiffusive and the second moment behaves as ~ A 2 NlnN for large N. 

When the discrete microscopic evolution equation (0) is approximated by a Fokker-Planck 
equation in the continuum limit, the elephant random walk can be shown to reproduce the 
dynamics of a Brownian particle in a time- varying harmonic oscillator potential Q]. This 
system is encountered in the study of anomalous diffusion that results from time- dependent 
Fokker-Planck drift coefficients jjJQ]. In this continuum approximation, an initially sharply 
peaked PDF of displacement asymptotically evolves into a Gaussian distribution with a 
time-dependent envelope. 

This Brief Report has two main contributions. First, we present a continuous time gen- 
eralization of the ERW by taking the time intervals between steps from a distribution for 
which the central limit theorem applies (Sect. This generalization leads to a more flexi- 
ble model as the waiting time between successive displacements (or events) in the ERW is 
no longer required to be constant. We will focus on an exponential waiting time distribution 
in a Poisson generalization that will allow us to obtain exact expressions for the moments 
of the distribution of displacement in terms of confluent hypergeometric functions (Sect. 
Illl|) . Second, we calculate the normalized third and fourth cumulants, or the coefficients 
of skewness and kurtosis, of the PDF of displacement in an attempt to discover new dy- 
namical transitions and to analyze the range of validity of the continuum Fokker-Planck 
approximation described above (Sect. ITV]) . 



II. THE CONTINUOUS TIME ERW 



We consider an elephant random walk in continuous time with a waiting time distribu- 
tion ip{t) that satisfies the following conditions: (i) This distribution represents mutually 
independent random waiting times, and (ii) It possesses finite first and second moments. 
For such a distribution, the probability density V(N,t) that N steps will take place in the 



ERW at time t is given by 



111 



v(N,t) = [m*] 



N 



i/j(t') dt' 



o 



(2) 



where the asterisk * is the Laplace convolution operator. Formally, the continuous time 
PDF of total displacement P(y, t) is therefore 

oo 

P(y,t) = J2'P(N,t)P N (y), (3) 

7V=0 

and its moments (y™) are 

oo 

(yr) = J2 V (N,t)(y%), (4) 

7V=0 

where (y 1 ^) are the moments of the N step PDF of displacement -P/v(|/) in discrete time. The 
central limit theorem guarantees that in the asymptotic limit the continuous time PDF 
P(y, t) becomes identical to the analogous distribution of the discrete time ERW having the 
same average waiting time. That is, if 1/A is the mean waiting time, then asymptotically 

P(y,t*N/\) ^P N (y) as N^oo. (5) 

Thus, all of the dynamical transitions observed in the discrete time ERW are also present 
in our continuous time generalization and occur at the same critical values. Furthermore, 
the sam e continue app—ions used in Q may als0 be appHed to this eontinuous t im e 
model with the same degree of success. 

From this well-known result, we are able to obtain exact closed form results for the 
moments of P(y, t) if we consider a waiting time distribution given by the exponential 

xfj(t) = \ e - xt . (6) 

The number density V(N,t) is now given by the Poisson distribution and the moments of 
P(y, t) are 

(yr> = e- At E^<yJ?>- (7) 

iV=0 



In the following section, we shall calculate the first four moments of P(y,t) in this Poisson 
ERW and show that they indeed approach the analogous discrete time moments as demanded 
by the central limit theorem JSJ). 



III. MOMENTS IN THE POISSON ERW 



To calculate the moments of P(y,t), we will need the exact expressions for the moments 
of the displacement in a discrete time ERW and substitute these values into (jZJ). The first 
two discrete time moments are given in as well as a general recursion formula that will 
allow us to calculate the higher moments. Quoting these results and solving for the third 
and fourth moments in an N step discrete time ERW yields 
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(11) 



Here, we have made extensive use of the rising factorial or Pochhammer symbols (a) at [i£IJ m 
anticipation of encountering hypergeometric sums as we solve for the moments of the PDF 
of displacement. 

To illustrate the evaluation of the continuous time moments, we shall calculate the mean 
total displacement explicitly. By making the appropriate substitution into (J7J) and consid- 
ering the action of the Xtd/d(Xt) operator, we are led to the expression 

(3Ae~ x \ ± d ^(Xt) N 



(yl) 



-Xt 



a 



d(Xt) 



E 

N=0 



N\ 



IN 



(12) 



Recognizing the power series representation of Kummer's confluent hypergeometric function 
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M(a, c; At) and performing the differentiation yield the following exact result: 

(yj) = f3AXte- xt M(a + l,2;Xt). 



(13) 



Similarly, we obtain the following expressions for the second, third, and fourth moments 
of the distribution of total displacement: 

A 2 At e~ xt r 



(y 2 t ) 



2a - 1 
/?A 3 At e~ xt 



2aM{2a + 1,2; Xt) - e xt 
6(a + l)M(3a + 1, 2; At) 



(14) 



(vt) 



2{2a - 1 

-2(a + 4)M(a+l,2;At) 

-3(a + l)AtM(a + 2,3;At) 

A 4 At e~ xt f 24a(2a 2 + 2a — 1) , n ^ 
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-4a(2a + 5)M(2a + l,2;At) 
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4a 2 + 12a - 1" 



(15) 



+ e 



3At 



(16) 



4a - 1 

Applying L'Hopital's rule as the parameter a is made to approach 1/2 allows us to find 
the scaling properties of these moments when the process is marginally super diffusive. As 
in the discrete time process, we find that the second, third, and fourth moments behave 
asymptotically as A 2 AtlnAt, 6vr- 1 / 2 /5A 3 (At) 3 / 2 In At, and 3A 4 (At In At) 2 , respectively At 
the critical point a = 1/4, a mild transition arising from the fourth moment is observed for 
small times but is obscured in the asymptotic limit. 

From the asymptotic series expansion of the confluent hypergeometric function we 
find that the leading terms in the asymptotic expansions of these continuous time moments 
are identical to that of the corresponding discrete time moments for At w N, as we asserted 
earlier. In Fig. ^ we confirm graphically that the dimensionless moments in a Poisson 
continuous time ERW approaches those in a discrete time ERW with the same mean waiting 
time. In these plots we have used several values of a to illustrate the behavior of these 
moments in several regimes of the ERW: superdiffusive (a = 1.0), correlated diffusive (a = 
0.499), no memory (a = 0.0), and anticorrelated diffusive (a = —0.5, —1.0). We find good 
agreement between corresponding moments as early as when t w 100/ A. 
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FIG. 1: (Color online) The moments of the total displacement in a Poisson continuous time elephant 
random walk (solid lines) asymptotically coincide with the moments of a discrete time ERW (dashed 
lines) with the same mean waiting time 1/A. Representative curves have the memory parameter 
a equal to 1.0 (•), 0.499 (o), 0.0 (▼), -0.5 (V), and -1.0 (■). Note that the odd moments vanish 
when the initial bias (3 is zero. 

IV. COEFFICIENTS OF SKEWNESS AND KURTOSIS 

In this section, we calculate the coefficients of skewness and kurtosis and examine their 
scaling properties in the long time limit. Solving for these quantities will yield relative 
measures of the skewness and kurtosis of the PDF in relation to its width and will allow 
us to identify dynamical transitions in these quantities. If we label the time dependent 
cumulants of the PDF of displacement as k™, the coefficients of skewness and kurtosis are 
defined as cf = /^/(k 2 ) 3 / 2 and c\ = k^/(k 2 ) 2 , respectively. 

By considering the leading order terms of the asymptotic expansions of these cumulants 
and taking the limiting value of c 3 and c\ as t/X approaches infinity, we find that these co- 
efficients vanish asymptotically when a < 1/2. This long time behavior under the condition 
of normal diffusive scaling is correctly reproduced by the continuum Gaussian approxima- 
tions discussed earlier and justifies its use in the diffusive regime. However, when a > 1/2 
the coefficients asymptotically approach constant values and the PDF of displacement is no 
longer approximately Gaussian. In this regime of superdiffusion the normalized coefficients 
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FIG. 2: (Color online) Asymptotic values of the coefficient of skewness (left) and kurtosis (right) 
in the superdiffusive regime of an elephant random walk. 

depend asymptotically on the memory parameters a and (3 according to 



cL = - / V^[r(2a ) r{ f ^ 

3(a + l) 2/? 2 (2a-l 



r(3a+l) [r(« + i)] 3 

,2 f (2a 2 + 2a - 1] 



(17) 



4 = -3 + 6[r(2a)] 2 |- 



(4a - l)r(4a) 
2/3 2 (2a 2 + a - 1) 2/5 2 (2a - 1) 



T(3a + l)r(a + l) T(2a)[r(a+ l)] 2 

" [T(a + !)]*/■ (18) 
These asymptotic values are plotted in Fig. 121 where we have labeled the critical contours 
where these coefficients vanish. 

We find that the skewness of the PDF has the opposite sign as the initial bias (3. Due 
to strong positive correlations in the superdiffusive regime, taking steps in the direction of 
the initial bias becomes even more favorable as the walk progresses. Thus, when the walk 
is initially biased in the positive direction, the most probable displacement becomes greater 
than the mean displacement at large times and the PDF becomes negatively skewed. 

There are two transition curves separating asymptotically platykurtic (c^ < 0) and lep- 
tokurtic (c^j > 0) PDFs shown in Fig. |21 The peak of the PDF of displacement becomes 
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relatively sharper than a similar Gaussian distribution as the dynamics of the walker be- 
comes ballistic when a approaches unity and (3 approaches ±1. Indeed, recent numerical 
calculations of the PDF of displacement in superdiffusive ERWs with bias (3 = 1 ^| in- 
dicate the presence of such deviations from a Gaussian distribution in the critical central 
region. The opposite behavior is observed as the initial bias (3 vanishes. We find that in 
the superdiffusive regime the memory effects result in a flatter distribution of displacement 
than the corresponding distribution of the memory-less Bernoulli walk Q|. 

The general non- vanishing skewness and kurtosis of the PDF and the additional memory- 
induced transitions in the superdiffusive regime that we have described above are not ob- 
served in the time-dependent Gaussian approximation obtained from a Fokker-Planck equa- 
tion. The origin of the breakdown of this approximation when a > 1/2 becomes clear if 
we consider a theorem by Pawula that states that a necessary condition for the jus- 
tified truncation of a Kramers-Moyal expansion of the microscopic evolution equation is 
that normalized cumulants beyond the second should vanish [l^. This condition is clearly 
not satisfied in superdiffusive elephant random walks. Although a Fokker-Planck approx- 
imation to the evolution equation will yield correct asymptotic mean and mean squared 
displacements, nontrivial information on the PDF of displacement that is related to higher 
non- vanishing cumulants is lost. 



V. CONCLUSION 



In this paper, we present a continuous time random walk model with complete memory 
of its history and calculate the exact moments of the distribution of displacement for the 
case of an exponentially decaying waiting time distribution. In this Poisson generalization, 
we show that the continuous time moments are asymptotically equal to the moments of the 
analogous discrete time process with the same mean waiting time as required by the central 
limit theorem. An analysis of the normalized third and fourth cumulants of the PDF of total 
displacement reveals new memory-induced transitions that appear in a parameter regime 
where this process exhibits anomalous diffusion (a > 1/2). The facts that these cumulants 
do not vanish and that these transitions exist show that Fokker-Planck approximations are 
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inadequate in the superdiffusive regime of an elephant random walk. 
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